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Abstract 

Locally relevant conditions, such as water stress in irrigated agricultural regions, should be considered when assessing the 
risk of crop allele introgression into wild populations following hybridization. Although research in cultivars has suggested 
that domestication traits may reduce fecundity under water stress as compared to wild-like phenotypes, this has not been 
investigated in crop-wild hybrids. In this study, we examine phenotypic selection acting on, as well as the genetic 
architecture of vegetative, reproductive, and physiological characteristics in an experimental population of sunflower crop- 
wild hybrids grown under wild-like low water conditions. Crop-derived petiole length and head diameter were favored in 
low and control water environments. The direction of selection differed between environments for leaf size and leaf 
pressure potential. Interestingly, the additive effect of the crop-derived allele was in the direction favored by selection for 
approximately half the QTL detected in the low water environment. Selection favoring crop-derived traits and alleles in the 
low water environment suggests that a subset of these alleles would be likely to spread into wild populations under water 
stress. Furthermore, differences in selection between environments support the view that risk assessments should be 
conducted under multiple locally relevant conditions. 



@ PLOS | one 



Citation: Owart BR, Corbi J, Burke JM, Dechaine JM (2014) Selection on Crop-Derived Traits and QTL in Sunflower {Helianthus annuus) Crop-Wild Hybrids under 
Water Stress. PLoS ONE 9(7): e102717. doi:10.1371/journal.pone.0102717 

Editor: William J. Etges, University of Arkansas, United States of America 

Received March 3, 2014; Accepted June 21, 2014; Published July 21, 2014 

Copyright: © 2014 Owart et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by the United States Department of Agriculture Biotechnology Risk Assessment Program competitive grant #2010-33522- 
21 668 from the United States Department of Agriculture - National Institute of Food and Agriculture. The funders had no role in study design, data collection and 
analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* Email: dechaine@cwu.edu 



Introduction 

Although gene flow from cultivated to wild populations has 
likely been occurring since the domestication of wild lineages, 
interest in the topic has increased with the commercialization of 
transgenic crops. Consequences of crop-wild hybridization may 
include the escape of engineered genes into wild populations [1], 
increased invasiveness of wild relatives of the cultivar [2] [3] , and 
the potential for crop-wild hybrids to outcompete native taxa [4]. 
Crop-wild hybridization must thus be considered when discussing 
the potential impacts of transgenic cultivars. The selective 
advantage of an allele is the best predictor of its establishment 
and spread in a new population [5] [6], and potentially 
advantageous transgenic alleles have been documented in 
Cucurbita pepo [7], Helianthus annuus [8], and Oyrza saliva 
[9]. These studies suggest that transgenes can increase fitness in 
crop-wild individuals following introgression, but the ubiquity of 
their effects on wild populations under diverse selective environ- 
ments is largely unknown. 

Although much attention has been given to the escape of 
engineered transgenes, crop-wild hybridization may also contrib- 
ute advantageous natural alleles to wild populations and can thus 
be used to model how selection acts on potential targets for genetic 
engineering. Hybridization between non-transgenic cultivars and 
wild relatives has resulted in the generation of at least seven types 



of agricultural weeds [10], and range expansion of wild 
populations following crop-derived allele introgression has been 
observed in Sorghum halapense, Rhododendron ponticum, and 
Manihot replans [11]. Studies in sunflower have suggested that 
crop-like traits (i.e., earlier flowering time and a larger primary 
inflorescence) increase reproductive output in crop-wild hybrids 
under a range of natural conditions [12] [13]. Crop-derived alleles 
have been known to persist within wild populations for at least five 
generations following hybridization in sunflower [14] and ten 
generations in wild radish [15], suggesting that some crop-derived 
alleles may persist in the wild, contrary to expectation. Indeed, 
crop-derived alleles may have also contributed to the evolution of 
weediness in wild sunflower [16]. Estimations of selection for crop- 
derived traits and alleles in wild environments may be applied 
toward understanding the consequences of transgene escape into 
the wild. 

The fitness effects of a trait or allele likely differ with 
environmental factors. For example, when exposed to increased 
interspecific competition and herbicide application, the relative 
fitness of sunflower crop-wild hybrids increased compared to their 
wild counterparts, suggesting that crop-like traits are more 
advantageous under certain conditions [17] [18]. Studies in 
sunflower also show that crop-like flowering (early) was advanta- 
geous in the absence of herbivory [12] [13], but wild-like flowering 
(later) was favored when pre-dispersal herbivory was considered 
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[13]. Increased pre- and post-dispersal predation has been 
observed for F, crop-wild hybrids relative to wild sunflower 
individuals [19] [20], suggesting that the presence of such 
herbivores could reduce hybrid fitness. The spread of otherwise 
advantageous crop-derived alleles may thus be mitigated by 
natural environmental factors. As such, fitness should be assessed 
over a range of locally relevant conditions. 

Water stress is an important environmental factor influencing 
the fitness of crop-derived traits and alleles that has not been 
studied in crop-wild hybrids. Under drought conditions, plants 
may be more susceptible to disease and consequently derive 
greater benefit from transgenic disease resistance [21]. Given that 
there will be heightened demand for transgenic drought tolerance 
technologies since irrigated agriculture will increasingly be 
operating under water scarcity [22], it is important for crop-wild 
hybridization risk assessments to identify crop-like traits conveying 
success under water stress in natural habitats. 

Characteristics that result in the efficient utilization of leaf 
surface area and corresponding leaf water potential at turgor loss 
are likely to be advantageous under arid conditions [23]. 
Cultivated sunflowers typically do not display these characteristics 
as selection for growth in water abundant agricultural environ- 
ments has resulted in a trade-off with drought tolerance [24]. In 
cultivated sunflower lineages, water stress reduces stem height, 
stem diameter, number of leaves, and leaf area [25] [26]. 
Decreased shoot and inflorescence mass have been observed in 
both cultivated and wild sunflowers when exposed to water stress, 
but cultivated individuals were more susceptible to drought in 
terms of wilting morphology [27]. Nevertheless, it is possible that a 
subset of crop-derived traits and alleles may confer an advantage 
in arid environments following introgression into wild populations. 

The application of quantitative trait locus (QTL) mapping to 
crop-wild systems can elucidate the fitness effects of cultivar- 
derived alleles in natural environments. Estimates of crop-derived 
allele additive effects among loci from QTL mapping, coupled 
with phenotypic selection analyses, improve our ability to predict if 
crop-derived alleles will introgress into wild populations. Crop- 
derived alleles with strong additive effects that affect traits 
conferring a fitness advantage under diverse wild-like environ- 
ments have a high probability for introgression. Comparisons of 
significant QTL may be made across populations and environ- 
ments, allowing for the identification of crop-derived alleles that 
are likely to be favored across a wide range of relevant natural 
conditions [28]. 

Here we assess selection on crop-like traits in sunflower crop- 
wild hybrids exposed to water stress. We then infer patterns of 
selection on QTL by predicting how selection may affect crop- 
derived alleles in the tested environments. Specifically, we use a 
low water treatment in comparison to a control water treatment to 
investigate: (i) the direction and magnitude of selection acting on 
crop-like traits in each environment; (ii) the genetic architecture of 
these characteristics; and (iii) the inferred fitness effects of crop- 
derived alleles under these conditions. 

Materials and Methods 

Study System 

Common sunflower (H. annuus L.) is one of five highest 
production oilseed cultigens worldwide [29] and a model for 
studying crop-wild gene flow. Its weedy, self-incompatible wild 
relative (H. annuus var. annuus) is native to North America and 
has a rangeF extending throughout that of cultivated sunflower 
[30] [3 1] with abundance in the central and western United States 
[32]. Centuries of artificial selection have resulted in morpholog- 



ical divergence. Differences between cultivated and wild sunflower 
are readily observed in the number and size of heads, branching, 
leaf shape and size, and seed (achene) size [33]. Cultivated and 
wild sunflower populations frequently hybridize [34] [35], and 
crop-derived alleles are known to persist in the wild [14]. 

Mapping Population and Experimental Design 

As previously described [33] [12] [36] recombinant inbred lines 
(RILs) were developed from a cross between an oilseed cultivar 
(cmsHA89, PI650572) and a wild H. annuus var. annuus 
individual (Ann 1238, PI659440) from Keith County, NE, USA. 
Seed for the present study came from the F 7 -F 9 generations of 146 
RILs from this population. 

Between 21-24 May 2012, eight seeds per experimental entry 
(RILs and parents) were germinated in a greenhouse at Central 
Washington University (CWU), Ellensburg, WA before trans- 
planting into the field. Seeds of all lines were either sown direcdy 
into biodegradable pots (Jiffy, 3"x3") with field-collected soil at a 
depth of 2.5 cm (91 RILs) or, in the case of lines with known poor 
or undetermined germination rates, subjected to a dormancy- 
breaking treatment prior to planting (55 RILs). The soil type was a 
mitta ashy silt loam with a high available water capacity [37]. The 
dormancy-breaking treatment involved starting seeds that had 
been sterilized with 3% hydrogen peroxide on filter paper in Petri 
dishes, treating them with 25 ppm Ethephon for 24 hours, and 
transplanting them into the greenhouse upon cotyledon emer- 
gence. 

Upon emergence of the first true leaves, 18-28 June 2012, 
individuals of the 58 RILs with at least four replicates per line 
were transplanted into the field at CWU using a randomized 
complete block design. Replicates of each RIL were evenly split 
among four blocks and fully randomized across 24 paired rows 
within each block. Spacing was 0.3 m between paired rows and 
0.4 m between plants within a row. The field was tilled before 
transplanting to simulate disturbed habitat that wild sunflower 
populations frequendy colonize. The field site was located from 
the nearest wild sunflower population by ca. 40 km, thereby 
ensuring isolation. This region, which is an important irrigated 
agricultural area, receives less than 5 cm of precipitation and has 
daytime relative humidity of ca. 45% during the June - September 
period when this experiment took place. 

To help ensure establishment, seedlings were watered twice 
weekly before beginning water treatments on 24 July 2012. At that 
time, two blocks each were assigned to control water (CW) and low 
water (LW) treatments, to simulate regional cultivated and wild 
sunflower water conditions, respectively. CW blocks were spray 
irrigated for two hours once a week, and LW blocks were not 
watered throughout the growing season. Irrigation ceased in the 
LW blocks pre-flowering to stress plants during reproductive 
growth phases. Hand weeding was performed to minimize the 
effect of asymmetric interspecific competition between treatments. 

Plant Characteristics 

Plants began flowering on 1 August 2012. Flowering day was 
recorded upon opening of the first disc floret on the primary 
inflorescence of the first flowering plant. At peak flowering, plant 
water status was collectively evaluated for each individual using 
leaf pressure potential and leaf water content. The largest leaf of 
each plant was collected pre-dawn on 9 September 2012 and used 
to measure leaf pressure potential in the laboratory using a 
standard pressure chamber (Soilmoisture Equipment Corp.). Leaf 
water content was measured on the same leaf as a proportion of 
dry to wet mass; each leaf was massed before and after drying for 
three days in an oven at 32°C. 
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At the start of senescence (sunflower developmental stage, R-8) 
[38], approximately 15 weeks of age, morphological characteristics 
were measured. Stem diameter was measured at the base of the 
individual and plant height was measured from the base to the 
capitulum along the primary stem. Petiole length, leaf (blade) 
length, and leaf width were recorded for the largest leaf on the 
plant. Head diameter was determined for the primary inflores- 
cence. Head total refers to the number of developing capitula 
exceeding 3 cm in diameter. Branch number was counted as the 
number of axillary shoots developed from the primary shoot. All 
seeds were collected and weighed for each individual. Total seed 
mass and the mean mass of ten seeds were used to indirecdy 
estimate total fecundity (i.e., total number of seeds produced). Leaf 
size was estimated by multiplying leaf length by width. 

Quantitative Genetic and Selection Analysis 

Only RILs with data from at least two replicates in each water 
treatment for all characteristics were included in analyses (except 
for the selection analyses), resulting in a total of 237 individuals 
from 32 RILs (see discussion). All analyses were conducted using R 
[39], except for QTL mapping. A restricted maximum likelihood 
(REML) approach was used to test for the fixed effects of 
treatment, block, and planting date, as well as the random effects 
of line (RIL) and line-by-treatment on each characteristic (lme4 
package) [40]. Models including only one random effect were 
compared separately (RLRsim package) [41]. Mixed model F- 
statistic significance scores were also produced (lmerTest package) 
[42]. All characteristics met model assumptions. Models were then 
used to calculate least-squares means (lsmeans package) [43] and 
variance components; broad-sense heritability was considered as a 
proportion of additive genetic (RIL) to total variance. REML 
models generated best linear unbiased predictors (BLUPs) used to 
estimate bivariate genotypic Pearson correlations within each 
treatment. 

Aster models for life-history traits were used to characterize 
phenotypic selection on characteristics within each water treat- 
ment (aster package) [44]. These models utilize maximum 
likelihood linear methods and improve upon previous least squares 
methods by modeling multiple components of fitness into a single 
variable, as well as by specifying a particular distribution for each 
component of fitness [45] [46] [47]. All characteristics were 
regressed with the combined fitness variable based on two 
components: survival to reproduction (Bernoulli distribution) and 
fecundity (truncated Poisson distribution). In order to account for 
selection on plants that died prior to reproduction, all individuals 
that survived to the initiation of watering treatments were included 
in selection analyses. Final sample sizes were 476 individuals from 
62 RILs in the CW treatment and 67 RILs in the LW treatment, 
respectively. Relative fitness was calculated as a proportion of 
mean population fecundity in each treatment [48]. Block was 
included as a model parameter. To test for significant effects of 
each trait, sub-models, each omitting one characteristic (df= 12, 
each), were compared against the full model (df= 13) using 
likelihood ratio tests. We also used aster to calculate 95% selection 
gradient (P) confidence intervals for each trait [49]. Variance 
inflation factors (VIF) [50] were used to assess whether any two 
characteristics violated assumptions of multicollinearity; however, 
no combination achieved a VIF greater than five and all 
characteristics were thus retained in the selection model. 

Crop- or wild-like designations were assigned to either positive 
or negative linear selection for vegetative and reproductive growth 
characteristics using previous trait assignments for these RILs [33] 
[13] and by comparing the mean cultivated parent trait value to 
the RILs. For example, when irrigated, the cultivar developed a 



July 2014 | Volume 9 | Issue 7 | e102717 



Selection in Sunflower Crop-Wild Hybrids under Water Stress 



Table 2. Character means and heritabilities. 





Mean (SD) 




Heritability (SE) 




Range (low, high) 






CW 


LW 


CW 


LW 


CW 


LW 


Fecundity 


365.56 (421.80) 


107.83 (127.10) 


0.24 (1.17E-4) 


7.33E-10 (2.15E-8) 


(0.33, 1238.00) 


(4.00, 309.80) 


Stem Diameter (mm) 


7.95 (4.74) 


5.34 (2.11) 


0.41 (9.02E-3) 


0.11 (9.55E-3) 


(2.84, 24.24) 


(2.57, 8.57) 


Plant Height (cm) 


56.98 (23.81) 


42.93 (15.30) 


0.64 (2.78E-3) 


0.37 (3.78E-3) 


(25.05, 126.83) 


(21.30, 70.24) 


Petiole Length (cm) 


5.19 (3.83) 


3.53 (1.97) 


0.24 (9.21 E-3) 


2.84E-10 (7.29E-7) 


(1.20, 14.28) 


(6.75, 88.20) 


Leaf Size (cm 2 ) 


82.1 3 (4.09) 


34.35 (31.21) 


0.73 (7.48E-4) 


7.81 E-10 (8.86E-8) 


(3.80, 20.87) 


(6.75, 88.20) 


Branch Number 


5.35 (4.59) 


2.22 (2.04) 


0.40 (1.04E-2) 


0.21 (2.38E-2) 


(0.33, 15.80) 


(0.00, 5.00) 


Head Diameter (mm) 


44.84 (21.05) 


37.56 (12.63) 


0.21 (2.48E-3) 


0.16 (2.88E-3) 


(13.88, 119.10) 


(18.50, 61.80) 


Head Total 


10.27 (9.04) 


5.28 (2.84) 


0.44 (5.64E-3) 


0.08 (6.76E-3) 


(2.40, 33.00) 


(1.50, 10.00) 


Days to Flower 


76.38 (12.72) 


77.29 (11.99) 


0.62 (6.38E-3) 


0.44 (5.41 E-3) 


(67.00, 1 1 2.40) 


(69.25, 1 1 1 .00) 


Leaf Pressure Potential (MPa) 


3.27 (1.37) 


4.41 (0.90) 


0.09 (1.16E-2) 


2.29E-11 (1.86E-7) 


(1.55, 6.55) 


(2.85, 6.00) 


Water Content 


20.91 (10.72) 


21.56 (7.55) 


1.70E-10 (1.06E-7) 


2.14E-10 (1.69E-7) 


(13.60, 34.30) 


(1 5.82, 34.40) 



Summary statistics for morphological, reproductive and physiological characteristics in sunflower recombinant inbred line (RIL) cultivar (cmsHA89) x wild (ann1238) 
hybrid populations grown in control (CW) and low water (LW) treatments. Least-squares means, standard deviations (SD), and heritability components were calculated 
using a restricted maximum likelihood analysis. Broad-sense heritability is the proportion of line (RIL) to total variance. 
doi:1 0.1 371 /journal.pone.01 0271 7.t002 



much larger head diameter than the RIL population (Table SI), 
supporting that a larger head is a crop-like trait. Few seeds of the 
wild parental population germinated, likely due to high dormancy, 
and data for it could not be obtained. Leaf pressure potential and 
water content could not be assigned a crop- or wild-like distinction 
because these characteristics were strongly influenced by water 
treatment (Table 1; see discussion). 

Genetic Mapping 

Total genomic DNA of the full set of 169 RILs was isolated 
from 200-400 mg fresh leaf tissue (142 RILs) using the DNeasy 
Plant Maxi Kit (Qjagen). For the 27 RILs that could not 
germinate, gDNA was direcdy extracted from 100 mg of seeds. 
Samples were initially ground in 2 mL tubes using a TissueLyser 
(Qiagen) before 1 mL of washing buffer was added to each sample 
[5 1] . Tubes were mixed on ice for 5 min and then centrifuged at 

Table 3. Phenotypic selection analyses. 



16,000 rpm for 10 min at 4°C. The supernatant was discarded 
and the pellet was used to extract gDNA following the DNeasy 
protocol provided by Qiagen. The quality of the DNA was 
evaluated by spectrophotometry using a NanoDrop ND-1000 
(Thermo Fisher Scientifc). To ensure that the quality of the DNA 
extracted from seeds and leaf led to equivalent results, 23 
randomly chosen RILs were genotyped using both types of tissues; 
no differences were observed between tissue types. 

An array targeting 384 single nucleotide polymorphisms (SNP) 
was developed from a subset of the SNPs present on the 10,640 
feature Illumina Infinium array described by Bachlava et al. [52]. 
SNPs were selected based on known polymorphism between the 
parents based on an initial screen with the larger array as well as 
known map positions [53]. Genotyping was based on the Illumina 
GoldenGate platform and the VeraCode technology (Illumina, 
San Diego, CA). Genotypes were manually called using Geno- 
meStudio V2011.1 (Illumina). A total of 353 SNPs were 





CW 






LW 








P 


Deviance 


P-value 


P 


Deviance 


P-value 


Stem Diameter 


(-3.08E-2, -2.72E-2) 


1 046.80 


<2.20e-16 


(-9.08E-2, -8.43E-2) 


2509.10 


<2.20e-16 


Plant Height 


(8.07E-4, 1.46E-3) 


45.82 


<2.20e-16 


(5.53E-3, 6.63E-3) 


473.78 


<2.20e-16 


Petiole Length 


(5.19E-3, 8.64E-3) 


66.95 


<2.20e-16 


(2.38E-2, 2.99E-2) 


232.64 


<2.20e-16 


Leaf Size 


(-6.36E-2, -4.28E-2) 


110.87 


<2.20e-16 


(2.13E-3, 2.54E-3) 


485.10 


<2.20e-16 


Branch Number 


(1.40E-2, 1.77E-2) 


290.06 


<2.20e-16 


(9.1 3 E-2, 9.85 E-2) 


2674.10 


<2.20e-16 


Head Diameter 


(1.14E-2, 1.18E-2) 


8246.50 


<2.20e-16 


(3.53E-2, 3.64E-2) 


17877.00 


<2.20e-16 


Head Total 


(4.01 E-2, 4.16E-2) 


10292.00 


<2.20e-16 


(-1.29E-3, 3.33E-3) 


0.75 


0.39 


Days to Flower 


(-1.31 E-3, -8.69E-4) 


92.84 


<2.20e-16 


(-9.97E-5, 2.08E-4) 


0.49 


0.49 


Leaf Pressure Potential 


(-1.12E-2, -6.85E-3) 


49.38 


2.1 1e-12 


(4.96E-2, 5.79E-2) 


852.04 


<2.20e-16 


Water Content 


(5.94E-3, 6.62E-3) 


1058.00 


<2.20e-16 


(1.21 E-3, 2.23E-3) 


43.57 


4.10e-11 


Results summary of aster model comparisons and 95% selection gradient ([)} confidence intervals for morphological, reproductive, and physiological characteristics in 
sunflower recombinant inbred line (RIL) cultivar (cmsHA89) x wild (ann1238) hybrids grown under control (CW) and low (LW) water treatments. Sub-models (d.f. = 12), 
each omitting one trait, were compared to the full model (d.f. = 13). Likelihood ratio test deviance and % z test P-values are shown. 



doi:1 0.1 371 /journal.pone.01 0271 7.W03 



PLOS ONE | www.plosone.org 



4 



July 2014 | Volume 9 | Issue 7 | e102717 



Selection in Sunflower Crop-Wild Hybrids under Water Stress 




PLOS ONE | www.plosone.org 



5 



July 2014 | Volume 9 | Issue 7 | e102717 



Selection in Sunflower Crop-Wild Hybrids under Water Stress 



Figure 1. Graphical representation of quantitative trait locus (QTL) composite interval mapping results from QTL Cartographer. 

Characters were mapped using cultivar (cmsHA89) x wild (ann1238) sunflower recombinant inbred lines grown in control water (blue) and low water 
(brown) treatments. Fill indicates the positive (closed) or negative (empty) additive effect of the crop-derived allele, and asterisks mark QTL with an 
additive effect in the direction favored by selection. 1-LOD (box) and 2-LOD (tails) thresholds are shown. 
doi:1 0.1 371 /journal.pone.01 0271 7.g001 



successfully scored and thus included in the analysis. Correlations 
between SNPs were estimated using the Pearson coefficient using 
the cor function [54] in R [39]. Visual examination of the 
correlation matrix as a heat map revealed that four SNPs 
(SFW00635, SFW08150, SFW00608, SFW05812) exhibited 
strong associations with distant positions in the genome. Given 
the history of duplication in the sunflower genome [55], this was 
presumably due to the mapping of different paralogs in different 
populations, as has been previously documented in sunflower [53]. 
Consequently, the relative positions of all SNPs for the 1 7 linkage 
groups were re-calculated using MAPMAKER/EXP 3.0 [56] [57] 
knowing a priori the assignment to each linkage group (LG). The 
Kosambi distance function was set for all linkage groups. The 
order of markers presumably belonging to a same LG was initially 
inspected using the "LOD" command. Based on this pairwise 
comparison matrix providing the distance (in cM) between SNPs 
and LOD scores, we made a core group of markers less than 5 cM 
apart and the relative order was evaluated by the "compare" 
command. The most likely order was kept and the remaining 
SNPs to be incorporated to the LG were sequentially added using 
the "try" command. At every step, the new order was evaluated 
using the "ripple" command. Once all SNPs were assigned to a 
LG, the map was enriched by adding 140 previously mapped 
simple sequence repeat (SSR) markers [58]. Each SSR was initially 
assigned to a LG using the "near" command. The relative position 
of the SSR within a LG was estimated using the "try" command 
and the final order was then checked using the "ripple" command. 

QTL Mapping 

QTL mapping was conducted using the SNP/SSR map on 
BLUPs for all characteristics using a composite interval mapping 
(CIM) protocol in QTL Cartographer [59]. Forward and 
backward regression settings were chosen (P — 0.05) with a walk 
speed of 2 cM. Significant logarithm of odds thresholds (LOD) 
were estimated with 1000 permutations tests per characteristic 
[60]. QTL Cartographer was used to calculate the additive effect 
of the cultivar allele and the percent variance explained by each 
QTL. The additive effect was standardized to one standard 
deviation. Multi-trait mapping to assess QTL x environment 
interactions was conducted in QTL Cartographer using afore- 
mentioned CIM settings [61]. 

Results 

Treatment Effects 

Significant genetic variation was observed (RIL effects) for all 
characteristics except water content and RILs differed in their 
response to the water treatments (RIL x treatment effects) for stem 
diameter, plant height, petiole length, leaf size, branch total, and 
head total (Table 1). For all characteristics, the degree of 
heritability was lower in the LW than the CW treatment (Table 2). 
Of note, the heritability of leaf pressure potential and water 
content in the LW treatment approached 0.00, suggesting little 
additive genetic control for these characteristics. 

With the exception of water content, watering treatment effects 
were significant (Tables 1). Individuals were generally smaller and 
had less reproductive output in the LW treatment, as well as 
displayed greater leaf pressure potential (Table 2). Despite 



significant treatment effects for most characteristics, the direction 
of significant among-trait correlations remained consistent across 
both watering regimes (Table S2). 

Selection Analyses 

Most characteristics indicated significant gradients for direc- 
tional selection in both treatments (Table 3). Many characteristics 
showed consistent selection across environments. In both treat- 
ments, greater fecundity (higher fitness) was associated with: 
smaller stem diameter and larger plant height, petiole length, and 
head diameter, as well as more branches and greater water 
content. In contrast, the direction of selection differed between 
treatments for leaf size and leaf pressure potential. In the CW 
treatment, plants had higher fecundity if they displayed smaller 
leaf size and lower leaf pressure potential; selection estimates were 
reversed for these characteristics in the LW treatment. A greater 
number of heads and early flowering day were favored in the CW 
treatment although selection estimates were non-significant for 
these characteristics in the LW treatment. 

Crop- or wild-like attributes may be assigned to each direction 
of selection for vegetative and reproductive growth characteristics. 
Three crop-like traits were favored by selection in each water 
treatment: head diameter and petiole length in both treatments, 
leaf size in LW and days to flower in CW. Due to strong influence 
of environmental variation, a crop- or wild-like attribute was not 
given to leaf pressure potential or water content. 

Genetic Map Construction 

The genetic map was composed of the expected 17 linkage 
groups and included a subset of 358 of the 384 targeted SNPs and 
a set of 140 previously described SSR markers [58] [12] [13], 
resulting in a total of 498 markers. The final map covered 
1968.2 cM (range =83.0-161.0 cM per linkage group), on par 
with previous studies in sunflower [62] [63], with an average 
intermarker interval of 3.9 cM (range =0.0-27.1 cM). 

QTL Mapping 

Despite the low sample size, 44 QTL affecting ten character- 
istics were detected across both treatments (Table 4; Fig. 1). 
Twenty QTL for ten characteristics explaining 10.5-55.3% of the 
variance and 24 QTL for nine characteristics explaining 10.5- 
37.1% of the variance were mapped in the CW and LW 
treatments, respectively. At least one QTL was detected in both 
treatments for most characteristics, except plant height and leaf 
size in the CW treatment, head total in the LW treatment, and leaf 
pressure potential in both treatments. QTL x environment 
interactions were not statistically significant for any characteristic. 
Although QTL mapping in the present study utilized a small 
number of RILs, the number of QTL detected for traits that have 
been previously mapped using the same RIL population were not 
considerably fewer than in prior studies [12] [13] (see discussion), 
and 40% of these were detected on the same linkage group as 
previously reported for the same traits (Fig. 1). 

Similar clustering of QTL on LG 1 7 for seven characteristics 
was observed in both treatments. The additive effect of the crop- 
derived allele did not always confer the more crop-like trait 
(Table 4). For example, the crop-derived allele at two QTL (LG 7, 
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j> | | 16) increased branch number, even though domesticated individ- 

I ■§ ^ 3 uals have been artificially selected to display limited branching 



£ " o 

S o co g phenotypes. In addition, cultivars had lower water content on 

£ o -S '5 average compared to hybrid individuals in LW treatment (Table 

S £ "§ a SI), despite detecting the crop-derived allele increasing water 

~s «^ content at one QTL (LG 1 5) in the LW treatment. For 

=L ti >, § characteristics affected by multiple QTL, only branch number 

5 and days to flower for the CW treatment, as well as fecundity in 

"g M ,1 „ the LW treatment, indicated consistent direction of additive effects 

_ :§ ■§ <= of the crop-derived allele at all QTL. Additive effects of the crop- 

B -o ° I derived allele were in the direction of phenotypic selection analysis 

js for 64% of QTL overall, although a greater proportion of crop- 
derived alleles were predicted for selection in CW (75%) vs. LW 
(55%). 



E S 
§ 2 

at ru 
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Whether crop- or wild-like, the majority of traits conferring a 
fitness advantage in both environments generally increased plant 
competitive ability regardless of water stress. Characteristics 

8. £ c < related to competitive ability play a particularly important role 

5 S. "g -d in stressful, resource-limited environments as empty niche space 

" S Si determines habitat availability [64]. Stress alters the dynamics of 

S 8 ~° I. biotic interactions [651 and accentuates the benefit of traits such as 

c rs o ™ early vegetative growth [66]. Some competitive traits, such as 
resource allocation to increased vertical growth in sunflower, that 
may provide an advantage for procuring establishment in 

g = 5 disturbed habitats against other colonizing plant species have 
been actively selected against in cultivars, where in an agricultural 

i gi^ m environment their proximity to other individuals is artificially 

ii _o a | manipulated [67]. Nevertheless, we found some evidence that 
> g 5--"= crop-derived alleles contribute to competitive growth, suggesting 
"5 .5 > c that they are likely to aid in the expansion of hybrid plants in wild 
5 — a "8 environments. 

2! § Although leaves in the LW treatment were comparatively 

g -o smaller, crop-like (i.e. larger) leaf size was favored by selection only 

!» E « in the LW environment. Reduction in leaf size is a common 

■- ^ a) _i response to water stress as it may be necessary to decrease resource 

S EgO l oss anc j achieve equilibrium for maximal potential growth in 

5 g -g £ relation to available resources [68]. Leaf size reduction is also 

g jd g <S important to increase water use efficiency [69] ; however, too great 

"Z ° ~ a reduction may limit the photosynthetic ability in Helianthus 

§ 3 c species [70]. In addition, smaller leaves with a reduced boundary 

1 J; 3 layer will not lower leaf temperature and thus will not create a 

c I* 2 8 water use efficiency benefit [71]. Greater petiole length may 

_2 | jg j influence the relative advantage of leaf size as longer petioles 

| c -a S promote leaves to twist away from direct sunlight [7 1] . It is difficult 

oi c » J to generalize our leaf size results, as water availability may not 

~S"i."5 always predict the same differential selection [72]. Cultivar- 

_ derived alleles led to reduced or increased leaf size in the LW 

g n > § treatment depending on the locus. Future studies should attempt 

S3 «t! to separate direct effects of leaf size on tolerance of water stress 

a; §■ .1 ™ i- from indirect effects due to its correlation with other traits, 

o | > o ? Selection also reversed across treatments for leaf pressure 

o §> § 2! £ potential. Positive selection for leaf pressure potential, as observed 

3 ~E a. S 5 m m e LW treatment, indicates that plants with a greater ability to 

o ^ -g £ move water through transpiration had greater fecundity. Unfor- 

§ "§ s tunately, through a combination of low heritability, greater 

~ m g~ o E influence of treatment than genotype, and no prior studies of this 

*j 1 < ^ ,2, characteristic in crop-wild hybrids, we are unable to assign a crop- 

™ | S £ like attribute for leaf pressure potential. Given that cultivated 

C 

rU 

5 U = ~ '— 



trj a x 



S | ~* 5 individuals display a greater degree of wilting relative to wild 
% ft = ^ '5 individuals when exposed to water stress [27], cultivars would 
presumably display greater susceptibility to water stress, in that 
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they would be less able than wild individuals to maintain leaf 
pressure potential as water availability is reduced. 

While the present study included fewer RILs, our QTL 
mapping results are comparable with previous studies that utilized 
a greater sample size in the same population [33] [12] [13]. Small 
sample size in combination with low heritability in mapping 
populations can result in an upward bias in the variance explained 
by each QTL [73]. The detection of a small number of QTL of 
large effect thus cannot be used as the rationale for the simple 
genetic architecture of traits. Another consequence of small sample 
size is the potential for the detection of false positive QTL [74], but 
the generation of permutated significant LOD values is an 
accepted method to account for this [75]. Importandy, the 
additive effect of the crop-derived allele in our study was in the 
same direction as predicted from the selection analysis for 
approximately half the QTL (excluding water content and 
fecundity QTL) detected in the LW treatment, suggesting that 
the repeated detection of selection favoring crop-derived alleles 
may be advantageous under a broader range of natural conditions 
than otherwise thought [14] [15]. 

For many QTL, although the additive effect of the crop-derived 
allele was in the direction of selection, it conferred a more wild-like 
phenotype. For example, although decreased leaf size is a wild-like 
trait, the crop-derived allele had a negative additive effect for three 
QTL affecting this trait across both treatments. This result is not 
uncommon in QTL studies [76]. In previous studies of the same 
sunflower population, mixed effects of crop-derived alleles were 
found for several morphological characters [33] [12] [13]. In 
lettuce (Lactua sativa x L. serriola), the allelic directions were as 
expected for QTL for major effect but mixed for QTL of 
intermediate effect [7] [78] [79]. Consequendy, failure to include 
genetic architecture methodology may present an incomplete 
consideration of crop allele introgression. 

One caveat to all QTL studies using RILs is that they are 
limited to allelic variation between the two parent individuals. 
Although the wild individual in our study was collected from a 
semi-arid site, its drought tolerance is unknown and may not be 
representative of wild sunflowers in general. Nevertheless, several 
studies have now found that crop-derived traits and QTL are 
predicted for selection in a wild environment [12] [13] [79], 
supporting the hypothesis that the ability for hybrids to colonize 
new ecological habitats may have arisen from the migration of 
crop-derived alleles into wild populations [80]. That some alleles 
confer an advantage in both cultivated and wild environments 
indicates that adaptability to cultivated environments is present in 
the natural variation of some wild populations [16]. Selection for 
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